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Abstract. Soft cyclo-synchrotron photons produced by 
thermal electrons moving in a magnetic field can be sub- 
sequently Comptonized to produce hard X-ray and 7-ray 
emission. An accurate determination of the amount of soft 
cyclo-synchrotron photons produced from such plasmas 
is clearly necessary. We present the complete solution to 
the cyclo-synchrotron problem. We calculate the emission 
produced by an isotropic distribution of thermal electrons 
moving in a magnetic field, for temperatures ranging from 
5 x 10 8 K to 3.2 x 10 10 K. The applicability of the standard 
synchrotron formula is evaluated, and we provide conve- 
nient fitting functions for the emission spectra within this 
range of temperatures. We discuss how the Comptonized 
spectra, for various temperatures and optical depths, are 
affected by the new results. 



Surprisingly, only some limiting results are available in 
this field. Previous results (Schwinger 1949; Rosner 1958; 
Oster 1960, 1961; Bekefi 1966; Pacholczyk 1970) have con- 
sidered extreme cyclotron or synchrotron emission where 
analytical formulae have been used. These calculations do 
not treat the transition from cyclotron to synchrotron ra- 
diation as the velocities of the electrons become quasirela- 
tivistic, and can drastically overestimate the total amount 
of emission produced unless we are able to solve the prob- 
lem exactly. Previous work in the mildly relativistic regime 
(Petrosian 1981; Takahara k Tsuruta 1982; Melia 1994) 
include approximations which are not assumed here (see 
Mahadevan, Narayan & Yi 1996). As far as we know a 
complete treatment of the problem has not be presented 
so far. 
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1. Introduction 

Radiation by quasirelativistic and fully relativistic elec- 
trons in magnetic fields is very common in astrophysics. 
Examples of sources include jets and lobes in radio galax- 
ies (e.g. Carilli et al. 1991), hot accretion flows onto neu- 
tron stars and black holes (e.g. Narayan, Yi & Mahade- 
van 1995), and relativistic fireballs in gamma-ray bursts 
(e.g. Paczyhski & Rhoads 1993; Meszaros, Laguna & 
Rees 1993). These are a few examples of where we en- 
counter thermal or non-thermal electrons radiating cyclo- 
synchrotron radiation in a magnetic field which is proba- 
bly of near-equipartition strength. Under certain circum- 
stances (e.g. Takahara & Tsuruta 1982; Narayan, Yi & 
Mahadevan 1995) this radiation can then be multiply 
Compton scattered to give rise to detailed high energy 
spectra from these objects. The ability to accurately cal- 
culate cyclo-synchrotron emission is clearly important to 
the study and understanding of these systems. 



We give an outline of the complete solution (see Mahade- 
van, Narayan & Yi 1996 for details). The power (en- 
ergy /time/steradian/frequency) emitted by an electron 
moving with a velocity parameter between /3 and /3 + d(3, 
in a frequency range du>, and at an observer angle 9, is 
given by (Schott 1912; Rosner 1958; Bekefi 1966), 
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8(y) is the Dirac delta function, J m (x) is the Bessel func- 
tion of order m, J' m {x) is the derivative of the Bessel func- 
tion, /3|| = [3 cos 6 p , and /3± = sin^p is the velocity param- 
eter parallel and perpendicular to the magnetic field, B, 
respectively. 7 = ^/l — /3 2 is the Lorentz factor. 

Each integer m in the summation in Eq.(l) corre- 
sponds to a harmonic. To obtain the emission in the cy- 
clotron limit (low /?), we expand the Bessel functions in 
Eq.(l) to lowest order and obtain (Schwinger 1949; Rosner 
1956; Bekefi 1966) 
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Fig. la shows an example with [3 = 0.3, which corresponds 
to this limit. The widths of each harmonic in Fig. la is due 
to artificially broadening the delta function in Eq.(l) to 
facilitate the numerics. 

In the opposite limit corresponding to the extreme syn- 
chrotron case (7^1), we are interested in large m's and 
Eq.(l) reduces to (Schwinger 1949; Bekefi 1966; Rybicki 
& Lightman 1979) 
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where x = u>/u> c , with u> c = (3/2)7 2 wj sin 6 p , F(x) is given 

by ' ' 



F(x) = x / K L (£)d(, 



(7) 



and Ks(^) is the modified Bessel function. Fig. Id shows 
an example with 7 = 20, which corresponds to this limit. 

These results are valid for electrons of a given veloc- 
ity. For a thermal distribution of electrons, we need to 
integrate the emission over a Maxwellian. This integral 
can be done in the synchrotron limit, for a fixed particle 
direction, sin^p, to give (Pacholczyk 1970), 
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where xm = 2u>/3u> 6 2 , 6 e = kT/m e c 2 , K^il/Oe) is the 
modified Bessel function, and I(X) is given by 
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Fig. 1. Emission produced by a single electron moving with 
(a) /3 = 0.3, (b) /3 = 0.6, (c) /3 = 0.9, (d) 7 = 20.0. Dashed 
curves represent the synchrotron approximation Eqn. (6). The 
x-axis is the log of the dimensionless frequency \ = w/w6, and 
the y-axis is r] x = rj^uib 

To obtain the exact emission spectrum for an isotropic 
thermal distribution of electrons, n(/3), for any tempera- 
ture, we evaluate L w = dE/dui given by 
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For a fixed [3 and u>, we numerically evaluate Eq.(10), as 
well as the sum over all the harmonics in Eq.(l). We repeat 
this calculation for various values of [3 and u> and tabu- 
late the results. The results can then be convolved with 
any isotropic velocity distribution to obtain the spectrum 
L w . In this paper we restrict ourselves to a relativistic 
Maxwellian n(/3) and present results for this particular 



3. Results 

3.1. Harmonic Nature of Synchrotron Emission 

Fig. 1 shows the total emission produced by a single elec- 
tron moving with different velocities. We plot rj x = dE/dx 
against x, where x = w/wj is the dimensionless frequency. 
In these units the vertical axis scales with the magnetic 
field, B, and the horizontal axis remains unchanged. The 
dashed curves are the synchrotron approximation Eq.(6). 
Even at high values of 7, the synchrotron approximation 
overestimates the emission at low frequencies x- This is 
due to discrete harmonic emission. 
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Fig. 2. Polar plots for a particle with 7 = 20, moving at an 
angle of cos 9 P = 0.3 to the magnetic field. The magnetic field 
is oriented along the positive y-axis. Different plots correspond 
to observing the emission at different frequencies (a) \ = 0.8, 
(b) x = 1-4, (c) x = 5.0, (d) x = 100.0. 

Fig. 2 shows polar plots of emission, at different x' s ; 
for a particle moving with a Lorentz factor 7 = 20.0 (cf. 
Fig. Id). At high values of x (Fig- 2d) we see the rel- 
ativistic beaming effect which is what we expect in the 
synchrotron limit. This is where the synchrotron approxi- 
mation and the numerical calculation agree. Discrete har- 
monic emission is dominant at low frequencies where we 
obtain "bumps" in the spectrum, and the agreement with 
the synchrotron approximation is poor (see Mahadevan, 
Narayan k Yi 1996). 



Fig. 3. Plots comparing the synchrotron approximation 
(dashed lines), Eqn. (8) (after averaging over particle direc- 
tions), with the exact numerical results (solid lines), for the 
emission from a thermal distribution of isotropic electrons at 
different temperatures, (a) T = 1.0 x 10 9 , (b) T = 2.0 x 10 9 , 
(c) T = 4.0 x 10 9 , (d) T = 32.0 x 10 9 . 

a, [3, 7, are adjustable parameters that are optimized so 
as to minimize the square of the deviation of M(xm) from 
the numerically calculated results. Table 1 shows the op- 
timized parameters for various temperatures. Fig. 4 com- 
pares the fitted functions (dashed lines), with the numer- 
ically calculated results (solid lines). The maximum error 
in the fits range from 10% at 10 9 A', to 0.1% at 3.2 x 10 10 A'. 



3.2. Integration Over a Thermal Distribution of Electrons 



We have computed emission spectra for isotropic thermal 
distributions of electrons with temperatures in the range, 
5 x 10 8 K < T < 3.2 x 10 10 K. The solid lines in Fig. 3 
represent the numerical calculations, and the dashed lines 
represent the synchrotron approximation, Eq.(8), after av- 
eraging over particle directions. It is evident that the syn- 
chrotron approximation for the emission overestimates the 
true emission at low temperatures, but is a good approx- 
imation to the true emission for T i> 10 10 K, as expected. 
We have obtained a set of fitting functions corresponding 
to each of the temperatures for which we have calculated 
the spectrum. The functions are of the form 
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with M(xm) replacing I(xm) in Eq.(8). With this sub- 
stitution Eq.(8) represents the total emission from an 
isotropic thermal distribution of electrons. In Eq.(ll), 
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Fig. 4. Plots comparing our fitting functions (dashed lines) 
with the exact numerical results (solid lines), for the emis- 
sion from a thermal distribution of isotropic electrons at dif- 
ferent temperatures. (a)T = 1.0 x 10 9 , (b) T = 2.0 x 10 9 , (c) 
T = 4.0 x 10 9 , (d) T = 32.0 x 10 9 . 
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4. Effects on the Comptonized Spectrum 

For certain temperatures and frequency ranges, Fig. 3 
shows how the synchrotron approximation overestimates 
the exact calculation. For plasmas where the optical depth 
r ^ 1 for all frequencies, this overestimate can drasti- 
cally effect the amount of Compton scattered photons in 
the X-rays and 7-rays. In plasmas where the emission is 
self absorbed for log(x) ^ 2, we are only interested in 
the cyclo-synchrotron emission for log(x) > 2 (Petrosian 
1981; Takahara & Tsuruta 1982; Narayan, Yi & Mahade- 
van 1995). Fig. 3 shows that the synchrotron approxima- 
tion is a fairly good approximation to the exact calcula- 
tion, in this regime, and the results presented do not effect 
the amount of Compton scattered photons. 

As an example we calculate the emission produced 
by accretion onto a 7.0 x 10 5 Mq black hole with M = 
lO~ 4 MEdd- We assume an isothermal sphere of radius 100 
Schwarzschild, with a uniform number density, magnetic 
field, and temperature. We calculate cyclo-synchrotron 
emission, Comptonization of the soft cyclo-synchrotron 
photons by the thermal electrons, bremsstrahlung emis- 
sion, and Comptonized bremsstrahlung emission. Fig. 5a 
shows the emission spectrum for a temperature of 10 9 K. 
The first peak corresponds to synchrotron emission, and 
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Fig. 5. Emission produced by accretion onto a 7.0 x 10 Mq 
black hole with M = 10" 4 M_B dd . (a) T = 1.0 x 10 9 , (b) 
T = 4.0 x 10 9 , (c) T = 8.0 x 10 9 , (d) T = 16.0 x 10 9 . 

the second to bremsstrahlung emission. Comptonization 
is not important at this temperature. In Figs. 5b, c, d, we 
see the Compton scattered emission, as the rising peak in 
the middle, as well as the Comptonized bremsstrahlung 
emission as a rise in the falling bremsstrahlung emission 
at the high frequency end of the spectrum. 



tion of thermal electrons moving in a magnetic field. For 
T ^ 3 x 10 10 K, standard synchrotron approximations 
overestimate the true emission. This discrepancy is im- 
portant when incorporating Comptonization for optically 
thin plasmas, since in some models (e.g. Narayan, Yi & 
Mahadevan 1995) it is the soft cyclo-synchrotron photons 
that are Comptonized by the thermal electrons. However, 
when self absorption is important, the Comptonized spec- 
tra are unchanged by the new results since the emission, at 
frequencies where the synchrotron approximation grossly 
overestimates the true calculation, is self absorbed. 
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Table 1. Optimal values of the parameters for different tem- 
peratures. 



T (K) 


I'(x) 




a 


P 


7 


5 x 10 8 


0.0431 


10.44 


16.61 


1 x 10 9 


1.121 


-10.65 


9.169 


2 x 10 9 


1.180 


-4.008 


1.559 


4 x 10 9 


1.045 


-0.1897 


0.0595 


8 x 10 9 


0.9774 


1.160 


0.2641 


1.6 x 10 10 


0.9768 


1.095 


0.8332 


3.2 x 10 10 


0.9788 


1.021 


1.031 
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5. Conclusion 



We have obtained the complete solution to cyclo- 
synchrotron emission produced by an isotropic distribu- 
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